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ABSTRACT 

In this talk, I will present a snapshot of key results regarding atmospheric 
neutrinos, as of the summer of 2002. I will review the evidence for neutrino 
oscillation, dominated by the large sample from Super-Kamiokande, but sup- 
ported by confirmation from MACRO and Soudan 2. I will then review some 
of the more detailed inquiries made using Super-Kamiokande data such as al- 
ternative scenarios, three-flavor oscillation, and tau neutrino appearance. 



1 Introduction 



Atmospheric neutrinos are produced by the decay of pions, kaons, and muons 
in cosmic ray showers. They have provided a fortuitous laboratory for neutrino 
oscillations, where the oscillation probability: 

P{va^ Vb) ^ sin^2esin^ (1) 



is probed over a wide range of L (10-10000 km) and E (0.5-100 GeV). 

Two normalizing principles are used to understand the data. First, cosmic 
ray showers consist mostly of pions, which decay to fi + v^, and the fj, decays 
to e + + j^e, resulting in a flux ratio {v^j^/v^) ^ 2. This ratio grows to larger 
values at high energy, as the time-dilated muon may strike the ground before 
decaying (removing Vf. from the beam), but at all relevant energies the ratio is 
predicted to about 5% accuracy. The second principle is that at equal zenith 
and nadir angles the flux of high energy neutrinos must be nearly identical. 
Again, this is predicted with an accuracy of a few percent. 

Before I review the experimental data, I should mention some develop- 
ments in the neutrino flux calculation 0) , which is a key input to more detailed 
understanding of the data. First, there is new data providing clarity regarding 
the flux of primary cosmic rays, which is the starting point for these calcula- 
tions. The BESS and AMS spectrometers have made new precision measure- 
ments 1,1) in the important range of 5-100 GeV/nucleon, with data in good 
agreement with each other whereas older measurements had fairly sizeable dis- 
agreements. 

Second, increased computational power and the large Super-K data sam- 
ple have motivated the transition from 1-dimensional calculations II) , where the 
Pt of the daughter particles is neglected and every neutrino is assumed to fol- 
low the trajectory of the primary proton, to 3-dimensional calculations S . In 
3-dimensions, one should account for pt as well as bending in the geomagnetic 
field. The first 3-D calculations showed a relatively surprising effect that turns 
out to play little if any role in oscillation studies. The effect is an enhanced 
neutrino flux at the horizon. Indeed, below ~ 1 GeV, there is fairly sharp 
peaking as a function of zenith angle, as shown in Fig. ||. The central dip in 
the 1-D curve is caused by geomagnetic cutoff; the central enhancement in the 
3-D curves overwhelms the cutoff. These effects are stronger with decreasing 
neutrino energy, but the lower the energy of the neutrino, the poorer is the cor- 
relation between the neutrino direction and that of the outgoing lepton. The 
net result is that the any structure at the horizon is washed out of the angular 
distributions of the lepton, as also shown in Fig. |l|. The predicted shapes are 
nearly the same for each calculation. There are some differences in the overall 
normalization between the models, but in neutrino oscillation studies, the flux 
normalization is usually allowed to float as a free parameter. 
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Figure 1: This figure compares 1-dimensional and 3-dimensional flux calcula- 
tions for the Super-Kamiokande site. The left panel compares the integrated 
V flux from 0.3 to 1 GeV as a function of zenith angle. The right panel shows 
the zenith distribution of /i-like data and Monte Carlo events in the Super-K 
analysis. The lower curves include zz-oscillation. 



2 Evidence for Neutrino Oscillations 

The evidence for neutrino oscillations §|) hinges on the two aforementioned 
normalizations of the atmospheric neutrino flux: that it should be composed 
of a well-predicted ratio of to v,, , and that it should be up-down symmetric 
at all zenith angles. To further cancel experimental systematic uncertainties, 
the flux ratio of z^^ to v^. is compared to a detailed Monte Carlo simulation by 
forming a double-ratio, R. For the sub-GeV sample of events {Ei, <, 1.5 GeV), 
the Super-Kamiokande analysis of 1489 days of data 0) shows: 

i? = i^^iZ^^il^ 0.688 ± 0.016 ± 0.050. (2) 

The up-down asymmetry is expected to be zero for the multi-GeV muon sample 
{El, >^ 1.5 GeV), but is measured to be: 

^ = ^^^^^^^^^^^ = -0.303 ± 0.030 ± 0.004, (3) 

whereas the asymmetry is consistent with zero for the multi-GeV electron sam- 
ple and the multi-GeV muon Monte Carlo. The up-down asymmetry deviates 
from expectation by more than lOcr, essentially independent of any Monte Carlo 
input. 



To estimate the parameters of neutrino oscillation, Super-K bins the 
single-ring Vf, and events in energy and zenith angle (effectively path-length 
L). Additional bins are also used: partially contained events where the 
muon exits the detector, multiple-ring events where the brightest ring is /i-like, 
or upward-going muons where the muon enters the detector from a neutrino 
interaction in the surrounding rock. Another category of multiple-ring events 
are selected where the brightest ring is e-like; these have an enhanced fraction 
of neutral current (~ 35%). These bins play a negligible role in the v^^ ^ 
fits but are valuable in restricting alternative models involving sterile neutrinos 
or neutrino decay. In total, 155 bins are formed, too many to easily display. 
Figure H shows 95 bins in 10 zenith angle distributions, where some of the 155 
bins have been combined. The data fits the prediction beautifully if a two- 
flavor v^J^ — oscillation is applied to the neutrino flux with sin^ W — \ and 
Ato^ = 3 X 10^^ eV^. The contours for this best fit are shown later, in Fig. |5[ 
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Figure 2: The zenith angle distribution for the atmospheric neutrinos used in 
the Super-Kamiokande oscillation analysis. 



In the past, the value of R measured by iron calorimeters such as NU- 
SExD and Frejus B seemed to show little evidence for neutrino oscillations, 
compared to the contemporaneous water Cherenkov detectors Kamiokande and 
1MB. Currently, the Soudan 2 fine-grained iron tracking calorimeter in Min- 
nesota, U.S.A., has recorded the largest data set 0) of contained atmospheric 
neutrino vertices using iron. They confirm the picture of atmospheric neu- 
trino oscillation, with some quite different systematics such as geomagnetic 
location, target nucleus, and reconstruction technique. For their analysis of 
5.9 kt-yrs of data, using a minimum momentum of 300 MeV/c, Soudan 2 re- 
ports R = 0.71 ± .09, in agreement with the results from water Cherenkov 
experiments. 

In addition, the Soudan 2 group exploits the good charged particle recog- 
nition of the detector to define a high resolution sample. The sample includes 
high energy quasi-elastics (single tracks), low energy quasi-elastics where the 
recoil proton from u -\- n —> I + p is identified, and high energy multi-prong 
events. In each case, the neutrino direction is well determined, with approx- 
imately 25° angular resolution (for z^^) and resolution in \og{L/E) of ~ 0.5. 
The data is binned by L/E and compared to Monte Carlo expectation for 
different oscillation parameters. Figure |^ shows the best agreement, found at 
1 X 10~^eV^ and sin^ 29 — 0.97. The apparent suppression of downward-going 
{\og{L/E) < 2.5) explains the high value of Am^ compared to the Super- 
K best fit value. When uncertainties are incorporated, the 90% confidence 
interval for this fit does allow Am^ down to a few times lO^'^eV^. 

The MACRO detector is too coarsely grained and lightweight (4.7 kt) to 
effectively study contained interactions. However, it has a solid angle accep- 
tance that rivals that of Super-Kamiokande, and is ideally suited to identify 
neutrino induced upward-going muons using time-of-flight. MACRO consists 
of a lower section filled with crushed rock absorber and an upper section that 
is hollow; both sections are instrumented with planes of streamer tube track- 
ing. Fast timing is measured by tanks of liquid scintillator that surround the 
detector, and the plane between upper and lower sections is also a timing plane 
of liquid scintillator. The simplest analysis is to measure the distortion in 
zenith angle flux for upward muons that pass completely through the detector; 
this is shown in the left panel of Fig. ^. These data demonstrate a best-fit to 
neutrino oscillation of 2.5 x lO^^eV^ and sin^ 29 — 1, in perfect agreement with 
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Figure 3: The L/E distribution for high rcsohition atmospheric neutrinos in 
the Soudan 2 analysis. 



the Super-K result. MACRO has also studied the approximate energy distri- 
bution of these events (in four bins of 14, 35, 90, and 150 GeV) by considering 
the amount of multiple scattering lil) . This analysis shows consistency with 
neutrino oscillations, with the greatest disappearance in the lowest energy 
bins. 

MACRO also measures the rate of lower energy 4 GeV) neutrino events 
that do not pass all of the way through the detector EH) . Internal upgoing 
events originate in the lower absorber but can be distinguished from downward 
stopping muons if they pass through two upper timing planes. In addition, 
there is a sample of internal downgoing events that start in the absorber and 
exit through the bottom. These have an indistinguishable topology from neu- 
trino induced upward stopping muons. Fortunately, both of these topologies 
are neutrino induced, with similar parent energy distributions. The results for 
these samples are shown in the right panels of Fig. ^. These data also favor 
neutrino mixing in the range of 10^'^ to lO^^eV^. 




Figure 4: The zenith angle distribution for atmospheric neutrinos of various 
categories in the MACRO oscillation analysis. 

Figure ^ shows the confidence intervals for the atmospheric neutrino anal- 
yses just described. The Super-Kamiokande, Soudan 2, and MACRO results 
are consistent with each other, with the Super-K result yielding the most sig- 
nificant and precise estimation of mass splitting and mixing angle. Despite the 
varied techniques of these three experiments, they do all rely upon atmospheric 
neutrinos, hence a further confirmation using a completely different neutrino 
source is desirable. The K2K experiment is the first long baseline experiment 
to probe neutrino oscillations with energy and fiight distance comparable to 
atmospheric neutrinos. K2K uses a 98% pure i/^ beam with mean energy 1.3 
GeV that travels 250 km from KEK to the Super-K detector. The preliminary 
results from about one half of the planned running have recently been made 
public 0) , and the contour from K2K is also overlayed in Fig. ^; the deepest 
minimum coincides exactly with the Super-Kamiokande best fit region. 
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Figure 5: Confidence intervals for the parameters of atmospheric neutrino oscil- 
lation, as determined by several different experiments. The two similar contours 
for K2K reflect different treatment of systematic errors in the analysis. 



3 Unlikely Alternatives 

An important question is whether the Super-K atmospheric neutrino observa- 
tions are only consistent with i/-oscillation or could alternatively be explained 
by neutrino decay lil) , or other exotic ideas such as violations of Lorentz 
invariance, flavor changing neutral currents etc. Table |l| summarizes the re- 
sults of a fit to the Super-K data, again binned in by energy and zenith angle 
for fully- and partially-contained events, upward going muons, and multi-ring 
events. The multi-ring events are subdivided into charged current z/^ and neu- 
tral current enhanced samples. The Monte Carlo sample used for i/^ ^ v^- 
includes tau lepton appearance. 

In Tabic |l] both the absolute and the differences in compared to the 
best fit of ^ Vt- oscillation are reported. In some cases, the functional form 



Table 1: Alternative fits to Super- Kamiokande atmospheric neutrino data. 
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of the transformation probability approaches that of the accepted hypothesis of 
pure ^ Vr- But in each case, there are typically several bins that constrain 
the alternative hypothesis. For example, matter effects suppress the oscillation 
probability for high energy /i-like bins (> 5 GeV x [Am^ in units of 10~'^eV^]). 
For another example, neutral current events should disappear with sterile neu- 
trinos or neutrino decay. By observing the differences in x^, one effectively 
ignores whichever bins are not sensitive to the differences in model. None of 
these pure alternatives fit the data very well; the decoherence model 13) is 
the most difficult to reject. Observation of an oscillatory minimum that is not 
washed out hy L/E resolution would greatly assist in discrimination. 

Clearly, one can generally include any exotic alternative as an admixture 
with "standard" neutrino oscillations. In fact, one is forced to consider both 
in the decay scenarios, where the neutrinos are required to have mass. For 
the important case of an admixture of sterile neutrino oscillation, Super-K 
has quoted limits on the allowable fraction of sterile neutrinos in a 4-ncutrino 
scheme li^. The v sterile mixture is found to be less than 20% (at 90% CL), 
with best fit to the Super-K data is at pure <-> Vt mixing. Figure |^ shows 
the allowed region a function of the admixture fraction. 

4 Three Flavor Oscillations 



Having established nearly maximal z/^ ^ Vr mixing as the explanation for 
the disappearance of atmospheric neutrinos, it is logical to search for a smaller 
component of mixing with v^- This is properly done in a 3-neutrino framework, 
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Figure 6: The allowed admixture of Vsteriie in a analysis of the Super- 
Kamiokande data. 



where a 3 x 3 unitary matrix mixes the mass eigenstates with the flavor eigen- 
states. Because the earth is full of electrons, but devoid of free muons or tauons, 
a component of electron neutrino flavor passing through the earth may be en- 
hanced by a matter effect resonance. The size of the resonant enhancement 
depends on neutrino energy and the pathlength through matter; therefore, this 
can potentially be observed by measuring a larger rate than expected in some 
narrow regions of zenith angle and energy. If one mass splitting dominates 
in scale over the other, as seems to be the case considering the likely results 
for solar neutrinos, the oscillation probabilities P(z/e <-> i/^), P(f^t <-> Vt), and 
P{i'r ^ ^e) may be expressed as functions of three parameters: two mixing 
angles (6*13, ^23) and the larger mass splitting Am^. 

The Super-Kamiokande group has analyzed the sample of electron neu- 
trino interactions with fine binning, but sees no deviation from the predicted 
rates. This null result is then turned into limits on the three mixing param- 
eters, as shown by contours in Fig. |^. The left panel shows that the small 
fraction of Vf. appearance allowed by the data is consistent with the limit on 
Vf. disappearance from the Chooz t2P reactor experiment. In fact, the data is 
consistent with no z/e mixing whatsoever. 
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Figure 7: Contours for the mixing parameters of three active neutrinos, in the 
single dominant Am^ scenario, as allowed by 1489 days of Super-Kamiokande 
data. 

5 Tau Appearance 

If Vt oscillation is taking place, tau leptons should appear in the at- 

mospheric neutrino flux. Most atmospheric do not have more than the 
« S.SGey of energy needed to produce a t lepton, however, in the large Super- 
K sample one does expect approximately 85 charged current Vr interactions. 
The events have many final state particles above Cherenkov threshold mak- 
ing exclusive reconstruction difficult, but a statistical analysis is possible. One 
preliminary analysis uses a neural network to select for characteristics of high 
pion production, based on event variables such as number of rings, number 
of ring seeds, number-of-decay electrons, particle classification of the brightest 
ring, and so on. 

For Am^ less than a few times 10~^eV^, it is expected that charged cur- 
rent tau appearance will only occur in the upward-going zenith bins. Therefore, 
the analysis is based on the up-down ratio of a sample selected for enhanced 
tau production, and the downward-going data may be used to check the per- 
formance of the neural network, as shown on the left of Fig. ||. This analysis 
has a 55% efficiency to save interactions in the fiducial volume while keep- 
ing 4.5% of the background. After cuts, a fit is done to the zenith angle 
distribution with components for both the background and the expected tau 
appearance signal. The data and fit results are shown on the right of Fig. H 




The best fit implies that the number of charged current Vr interactions was: 
99 ± 39stat ± 13Am2]!!5'5(3,y) ■ "^he first systematic error comes from varying the 
assumed Am^ between 1 and 5 x lQ~^eV'^ and the second comes from allowing 
three flavor oscillation at the CHOOZ limit. Although this is only a 2.3ct result, 
it shows that the Super-K data set is consistent with the expected rate of tau 
appearance. 
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Figure 8: The left panels show the performance of the neural network iden- 
tification. The right panel shows the zenith distribution of events selected as 
r-like; the filled region shows the relative proportion of t appearance that must 
be added to the atmospheric neutrino Monte Carlo to agree with the data. 



6 Super-Kamiokande Accident and Reconstruction 

On Nov 12, 2001 at 11:01 JST, a single photomultiplier on the bottom of the 
Super-Kamiokande tank imploded. At that time, the tank was 2/3 full of 
water; it was refilling after PMT replacement work in the preceding summer. 
The imploding tube caused a chain reaction that destroyed 6777 inner PMTs 
and 1100 outer PMTs. Subsequently, the collaboration made several studies 
to determine the mechanism of the chain reaction. A corrective measure was 
devised: the PMTs should be encased in a secondary shell consisting of a 



transparent acrylic hemisphere over the photocathode and a fiberglass shell 
surrounding the rest of the tube. This enclosure is not intended to withstand 
static water pressure; instead, several small holes are placed in the shell. Then, 
if a PMT were to implode, water rushes into the resulting vacancy more slowly 
and a high pressure shock wave is not developed. This has been tested several 
times in situ up to the full level of the Super-K tank. The tests show that if 
the center tube of a grid of 3 x 3 tubes is imploded, the damage is contained 
within the protective vessel, and the eight surrounding tubes do not implode - 
even if they themselves are unprotected. 

It will take several years to fabricate replacement photomultipliers for the 
inner detector. In the meantime, the detector is being recommissioned using 
the remaining tubes plus some spares to populate every other lattice point 
of the original photomultiplier tube array. This will be sufficient for nearly 
all research to continue, particularly the K2K long baseline experiment. The 
detector should be operational by the end of 2002. 
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